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ABSTRACT 

Using the Rossi X-ray Timing Explorer and the Nordic Optical Telescope we have ob- 
tained the highest ever quality X-ray/ white-light high-speed photometry of XB 1916— 
053. We refine the X-ray period (Px) to 3000.6 ± 0.2 s via a restricted cycle counting 
approach. Using our complete optical lightcurve, we have extended the optical period 
(Popt) ephemeris by another 4 years, providing further evidence for its stability, al- 
though a slightly longer period of 3027.555 ± 0.002s now provides a marginally better 
fit. Moreover, modulations at both Px and P opt are present in the optical data, with 
the former dominating the nightly lightcurves (i.e. a few cycles of data). We have also 
attempted to determine the "beat" period, as seen in the repeating evolution of the 
X-ray dip structure, and the variation in primary dip phase. We find that a quasi- 
period of 4.74 ± 0.05 d provides the best fit to the data, even then requiring phase 
shifts between cycles, with the expected 3.90 d "beat" of Px and P op t appearing to 
be less likely. Finally, considering the nature of each of these temporal phenomena, 
we outline possible models, which could explain all of the observed behaviour of this 
enigmatic source, focusing on which of Px or P opt is the binary period. 
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1 INTRODUCTION 

The low-mass X-ray binary (LMXB) XB 1916-053 was 
the first of the 10 dipping sources to be discovered. The 
recurrence period of the dip ~3000 s (Walter et al. 1982; 
White & Swank 1982) also provided one of the first bi- 
nary period determinations for a LMXB, and places it 
as the shortest period dipper. It also exhibits type-I X- 
ray bursts confirming that the compact object is a neu- 
tron star. A low- mass (~0.1 M Q ) donor is then required 
to fit into the system, which is most likely evolved, par- 
tially degenerate and hydrogen depleted (Nelson, Rap- 
paport & Joss 1986). The lack of eclipses by this com- 
panion, for a representative radius ~O.li?0, yields an 
inclination in the range 60° to 79° (Smale et al. 1988). 
The generally accepted model for the origin of X- 
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ray dips, suggests that they are due to the obscuration 
of the central source by vertically extended material at 
or close to the outer edge of the accretion disc (Par- 
mar & White 1988; Armitage & Livio 1998). In recent 
years Church and co-workers have studied a selection 
of dippers, utilising the superior spectral capabilities 
of ASCA and BeppoSAX, to probe in detail the spec- 
tral evolution through the dips (Church & Balucinska- 
Church 1993; Church & Balucinska-Church 1995). They 
find that a partial covering absorption model repro- 
duces the observed spectra well, where the fractional 
contribution from each of two components (point-like 
blackbody and more spatially extended cut-off power 
law) changes independently For XB 1916-053 a physi- 
cal model is implied consisting of the point-like black- 
body emission region (presumably the neutron star sur- 
face and/or adjacent boundary layer), plus the cut-off 
power law component, originating in an extended accre- 
tion disc corona (ADC, responsible for Comptonization, 
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Church et al. 1997). In the non-dip spectra the black- 
body alone contributes ~30% to the total flux (Church 
et al. 1997; Church et al. 1998). During dips the black- 
body component is very rapidly and totally obscured, 
whilst the power-law is more gradually attenuated, as 
the partial-covering fraction increases (even to complete 
coverage in the deepest dips seen). From measurement 
of the ingress/egress times a size for the radial extent 
of the ADC is found of 20% of the tidal disc radius. 
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Figure 1: Determinations of the "orbital" period of 
XB 1916-053. The X-ray period (-3000s) measure- 
ments from a long-line of satellite missions at the top 
and the three previous optical results at the bottom (l. 

White & Holt 1982, 2. Walter et al. 1982, 3. Smale et al. 1989, 4. 
Yoshida et al. 1995, 5. Church et al. 1997, 6. Church et al. 1998, 
7. Morley et al. 1999, 8. this paper, 9. Callanan, Grindlay & Cool 
1995, 10. Grindlay et al. 1988, 11. Grindlay 1989, 12. Schmidtke 
1988). 



1.1 Multi-periodicities of XB 1916 053 

Following their discovery of an optical counterpart to 
XB 1916-053 Grindlay (1988) searched for and found 
an optical modulation close to the 3000 s dip recur- 
rence period. Subsequently, this source has been the 
object of many timing studies. It was soon apparent 
that the periods derived from X-ray and optical obser- 
vations were inconsistent, differing by —1%. The results 



of these determinations are summarised in Figj^. For 
the X-ray dips, the earliest method adopted (by Wal- 
ter et al. 1982; White & Swank 1982), was to consider 
the timings of a series of (near-consecutive) dips from 
a single satellite pointing, and calculate the best fit re- 
currence time. Unfortunately the short observing base- 
line leads to large uncertainties in the period. More re- 
cently, the approach generally taken has been to epoch- 
fold the lightcurves on a grid of trial periods, and com- 
pute chi-squares to identify the period yielding the most 
distinct modulation. However, given the notable varia- 
tion in the morphology of the dip structure over time 
(see e.g. Smale et al. 1988; Smale et al. 1992; Yoshida 
et al. 1995; Church et al. 1997), the precision of this 
method is also uncertain, and the large errors quoted 
may still be underestimates. In their optical studies, 
Schmidtke (1988) and Grindlay (1988) used a straight- 
forward Fourier transform of their long (and irregularly) 
sampled lightcurves, whilst Callanan, Grindlay & Cool 
(1995) refined this by way of cycle-counting using the 
minima in the folded optical lightcurves to define the 
ephemerides. 

Not only are two close periodicities present, there 
is evidence for two longer periodic phenomena. In fact, 
it is assumed that physically one of the two ~ 3000s 
periods arises from binary motion, whilst the other is 
due to a beat with another longer period effect, which 
must be close to 4 days. Indeed changes in the lightcurve 
shapes in both bands and in particular the phases of 
dip features has been seen on this timescale (Grindlay 
1989; Smale et al. 1989; Yoshida et al. 1995). Note, for 
brevity, we shall henceforth refer to any ~4 d period as 
the "beat" period of Px and P op t , whilst reserving the 
notation Pb ea t for the quantity calculated as {PZ — 
P ~J) _1 , and using P' hf , at for any other period close to 
this, i.e. not assuming that they are the same. Lastly, an 
analysis of Vela 5B data by Pricdhorsky, Terrell & Holt 
(1983) yielded a possible longer-term variation in X-ray 
flux of 199 days. But with a false alarm probability of 
10-20% (Smale & Lochner 1992), this detection was by 
no means firm. 

A number of models have been proposed to explain 
the multi-periodicities, as we shall discuss in section [t], 
but no clear consensus has been achieved given the lim- 
ited constraints imposed by the observations to date. It 
was clear that further complementary X-ray and opti- 
cal observations were vital in order to solve the mystery 
of the multiple periodicities and constrain the range 
of proposed models to explain them. Hence, we ob- 
tained time with both the Rossi X-ray Timing Explorer 
(RXTE) and the Nordic Optical Telescope (NOT) to 
undertake a simultaneous X-ray/optical campaign. It is 
the results of these observations that we present here. 
In section |^ we summarise the observations and data re- 
duction performed, moving on to the temporal analysis 
in § H (X-ray period) , § || ( "beat" period) , § [| (longer- 
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term variations) and § || (optical periods). We discuss 
our results in section [?]. 



OBSERVATIONS AND DATA 
REDUCTION 



Table 1: X-ray and optical observations of XB 1916-053 



Instrument / 


Date (UT) 


Time 


span 


Band 




(h) 


(binary 








orbits) 


RXTE/PCA 


23.06.98 


0.64 


0.76 


X-rav f2-60 keV) 

J\. lay 1 Z; \J\J rvC V / 


24 06 98 


0.60 


0.72 




25.06.98 


0.60 


0.72 




26.06.98 


0.59 


0.70 




27.06.98 


0.27 


0.33 




17.07.98 


5.89 


7.07 




18.07.98 


8.95 


10.74 




1 9 07 98 


11.31 


1 3 57 




20.07.98 


9.65 


11.58 




21.07.98 


5.72 


6.87 




23a.07.98 


4.29 


5.15 




23b.07.98 


6.58 


7.89 




24.07.98 


4.04 


4.85 




25.07.98 


4.05 


4.86 




26.07.98 


3.82 


4.59 




01.08.98 


3.61 


4.33 




10.08.98 


2.85 


3.42 




14.09.98 


1.28 


1.54 




16a.09.98 


0.23 


0.28 




16b.09.98 


0.94 


1.13 


Nordic Optical 


25.06.98 


2.29 


2.75 


Telescope/CCD 


25/26.06.98 


1.94 


2.33 


white light 


26/27.06.98 


1.94 


2.33 




27/28.06.98 


1.94 


2.33 




29.06.98 


2.11 


2.53 


RXTE/ASM 


24.02.96 


1355d 




X-ray (2-10 keV) 


-09.11.99 







XB 1916-053 and its optical counterpart were ob- 
served partly simultaneously using RXTE and NOT, 
La Pahna in 1998 June. Unfortunately, the coverage 
available by the satellite was very restricted, and fur- 
ther extensive observations were obtained in the period 
1998 July 17-September 16. In addition we have consid- 
ered the long-term flux behaviour as measured by the 
All-Sky Monitor (ASM) on-board RXTE. The log of all 
these observations is given in table [l], whilst their rela- 
tive distribution and overlap is shown in Fig.|] (together 
with a plot of the the long-term ASM lightcurve) . 



2.1 RXTE/PCA 

For our short-term variability analysis we made use 
only of data taken with the Proportional Counter Array 
(PCA) instrument (Bradt, Rothschild & Swank 1993) 
on RXTE. This consists of five nearly identical Propor- 
tional Counter Units (PCUs) sensitive to X-rays with 
energy range between 2-60 keV and a total effective 
area of ~ 6500 cm 2 . The combination of large area and 
an energy range covering the peak of the X-ray spec- 
trum of XB 1916-053, means it provides good signal-to- 
noise lightcurves. Since the source is an X-ray burster 
the PCA data were obtained in two triggered high time- 
resolution burst modes, so that X-ray bursts could also 
be studied in detail (the subject of a future paper), as 
well as in both the standard modes. However, for our 
study of the periodic temporal phenomena presented 
here, we excluded the bursts. The 'standard2', 16s time 
resolution data were extracted using XSELECT, and then 
background-subtracted using the Epoch 3 bright source 
model. Given the average count-rate in the PCA of 
100-300 cts s _1 , this correction should be perfectly ad- 
equate. Lastly, the timing information was transformed 
into HJD, to match up with that of the optical data. 



2.2 NOT 

The optical counterpart to XB 1916-053 was observed 
using the Nordic Optical Telescope (NOT), La Palma, 
during 5 consecutive nights starting 25 June 1998. Each 
night the observations covered at least two orbital peri- 
ods (i.e. 2 hours minimum per night). The observations 
were obtained using ALFOSC, which is a multi-purpose 
medium/low resolution spectrograph and imager. AL- 
FOSC was equipped with a 2kx2k Loral-Lesser thinned 
AR-coated CCD detector, which was operated in a sub- 
windowed mode enabling a readout time of the order of 
3 s. This, coupled with exposure times of 10 s then yields 
a true time resolution of ~ 13 s. All the observations 
were made in white light for maximum signal to noise 
and best possible time resolution. Note the Loral-Lesser 
chip has an almost flat high (75-85% ) QE over most 
of the optical range. All the images were then bias sub- 
tracted and flatfielded in the usual manner. As the tar- 
get is located just ~2"away from a 2 mag brighter star 
good seeing is vital for the photometry. During our run 
the seeing was generally better than 1.0", which allowed 
us to use the IDL implementation of the DAOPHOT 
routines to reliably extract the light curves of the tar- 
get and the 2 comparison stars within the same frame. 
Differential photometry was then performed to obtain 
the final lightcurves. 



2.3 RXTE/ASM 

Unsurpassed long-term monitoring of the flux from all 
the bright X-ray sources is available from the All-Sky 
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Figure 2: Temporal overview of the various datasets used in this paper. Upper panel: the 10 d-averaged lightcurvc from 
the RXTE/ ASM with times of the pointed observations marked with symbols: solid triangles- X-ray (RXTE/PCA), 
open diamond- optical (NOT). The vertical arrows indicate the predicted times of maximum X-ray flux, according 
to the 83 day (quasi?) periodic modulation found in the long-term ASM X-ray lightcurve. Inset: detail of the X-ray 
lightcurves from all the PCA observations (points), together with the ASM monitoring (histogram). Lower panel: 
detail of the simultaneous observations in June 1998 with the optical lightcurvc shown uppermost (open diamonds) 
and the X-ray (PCA) below (solid triangles). 



Monitor (ASM) of RXTE . This consists of three wide- 
angle shadow cameras equipped with position sensitive 
Xenon proportional counters with a total collecting area 
of 90 cm 2 . They work in the 2-10 keV energy range, 
providing roughly a timing resolution of 90 minutes be- 
tween dwells (i.e. 80% of the sky is scanned every 90 
minutes on average), and are sensitive to sources as faint 
as 30 mCrab. 

We extracted the complete lightcurve of XB WIG- 
OSS (spanning 1355 days) from the HEASARC database 
(the definitive flux values for each dwell), which contains 
both total flux and that in three energy bands. The av- 
erage count rate is only about 1.2 s _1 (16 mCrab), hence 
rcbinning is required to obtain adequate signal-to-noisc. 
We derived both 1 day and 10 day averages, using the 
same screening as that performed by the MIT ASM 
team, but also a requirement that there were at least 
10 (100) dwells per bin. This latter criterion ensured 
all data points had approximately the same estimated 



uncertainty, and generally excluded poor data. How- 
ever, we also found it necessary to exclude (by hand) 
the block of data from December 1998 -January 1999, 
when the viewing of XB 1916-053 led to particularly 
poor data. The effect of restricted viewing was also ev- 
ident in a number of other galactic centre X-ray bina- 
ries; each exhibited significant and variable increased 
flux levels during these months. 



3 REFINING THE X-RAY PERIOD 

Given the well-documented morphology changes seen in 
the X-ray lightcurves of XB 1916-053 we firstly decided 
to examine the dip behaviour in our own RXTE/PCA 
data. However, due to the short time (< 2700 s) on 
source per satellite orbit we had to use the combined 
lightcurve from a few such segments and phase fold, to 
provide full phase coverage of the dip period. Hence, we 
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Figure 3: The PCA X-ray lightcurves folded on the refined Px = 3000.6s, with 40 phase bins. Each panel shows the 
data from a single pointing (which occurred at most daily). In two cases, the pointing has been first split into two 
subsets, as we noted dramatic changes in morphology during that time. For clarity two full cycles are shown, and 
different scalings of the count rate axes have been employed. Note the clearly evolving morphology, which returns 
to approximately the same state every 4-5 days. The lightcurves marked by * were used to define cphcmcridcs for 
cycle counting (see text). 



first assumed a representative X-ray period of 3000.0 s 
(see Fig. [l]), and folded the data from each day's ob- 
servations on this period. This procedure produced in 
most cases a good signal-to-noisc average lightcurve (af- 
ter phase binning), with complete phase coverage. The 
final set of resulting folded lightcurves (following the 
period refinement) are presented in Fig.^. 

The striking cyclical variation in morphology is im- 
mediately apparent. It is clear given the highly non- 
sinusoidal nature of the modulation that a Fourier- 
based period search is not appropriate. Furthermore, 
even the use of the phase dispersion minimisation 
(PDM, Stellingwerf 1978) pcriodogram may not be 
ideal. However, to make progress we did perform a PDM 
analysis and the results are shown in Fig.^. This yields 
Px = 3001.0 ± 3.0 s, where the precision is limited by 
the large morphological changes and the scatter they 
produce. 

To refine this result, we attempted for the first time 
a cycle-counting analysis of the dip timings. Firstly, 



we took the period found using the PDM, and then 
folded and phase-binned each day's data on this, with 
an arbitrary phase zero (To), to derive an average, good 
signal-to-noisc profile for the dip. Even by inspection it 
is clear that these phases vary systematically with the 
periodic overall shape changes. Hence, we decided to 
derive ephemerides for the principal dip only for folds 
which appear to be at least at roughly the same phase 
in the "beat" cycle (i.e. whenever we have a single nar- 
row dip morphology), as also indicated (by *) in Fig. |[ 
Next, a Gaussian profile was fitted to these dips to es- 
timate their centre phase (4> d i P ), with a generous error 
given by the FWHM of this fit (rather than the statisti- 
cal error on the centre value). Thence, the time for the 
dip minimum is given by T dlp . n = T^ n + (j> dlp x P dlp , 
where T^ n = T + int[(T midpt - T )/P dlp ] x P dip . in or- 
der to give a time close to the mid-point time (T midpt ) 
of the nth dataset. We then have a set of 6 ephemerides 
spanning the full 3 months, which should all agree with 
a single recurrence period. For each cycle and trial pe- 
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Figure 5: For the best fit period of 3000.6 s, the observed (O) - computed (C) times of dip minimum in out 
RXTE/PCA data (for narrow dips only) are shown, in terms of cycles. 



riod, a time of minimum is predicted ( O) which is then 
compared to the observed time (O), and a simple x 2 
minimisation for the O-C values versus trial period then 
gives the best fit period of 3000.6 ± 0.2 s. The plot of 
X 2 versus Ptriai (Fig. [|) does show that the data can 
be fit by a number of periodicities, but only this result 
is consistent with the foregoing PDM analysis (not to 
mention the published results, see FigM). The O-C plot 
for the best fit period is shown in Figlq, which demon- 



strates the very small scatter of the dip phase for this 
period, over the 3 month time span! 



4 INVESTIGATING THE "BEAT" PERIOD 

We have previously noted the evolving X-ray lightcurve 
morphologies, whose timescalc of repetition is approxi- 
mately that of the "beat" period. An estimate of such 
a period is important as we believe the changes to be 
directly related to a precessing disc or whatever pro- 



On the multi-periodicities in the X-ray dipper XB 1916-053 7 



vides the longer period clock in XB 1916-053. Clearly, 
the energetics and extent of the disc/accretion stream 
interaction region (which provides the obscuring mate- 
rial for the dips) is varying, even significantly within a 
single visit (June 23 and September 16). In principle 
detailed modelling of these lightcurves can tell us much 
about the time- varying geometrical structure of the disc 
and how the stream impacts it. However, in this paper 
we will simply concentrate on the basic timing aspects 
and draw qualitative conclusions (see § 0). 

With our significantly improved measurement of 
Px we may now calculate Pbeat = (Px 1 — -P pt) _1 to 
much higher precision. Adopting the published value of 
P opt = 3027.551 ± 0.004s yields P beat = 3.90 ±0.03. But 
is this period consistent with the repeating changes in 
lightcurvc morphology? In figure [| we have plotted all 
of the folded X-ray lightcurves, vertically offset accord- 
ing to their phase assuming a 3.90 d period. There are 
at least two cases where from one "beat" cycle to the 
next, totally different dip morphologies would appear 
at the same "beat" phase, and hence we exclude this 
period. 

We will now attempt to determine a "beat" period, 
which is in agreement with the morphology changes. 
However, considering when the narrow dips or any other 
identifiable forms repeat may only crudely constrain the 
period to between about 3.5 and 5 days. We need an- 
other more quantitative method to indicate candidate 
periods. In principle there arc a number of ways to 
do this, namely: area of dip (i.e. below non-dip level), 
length of a dip, depth of dip, and phase of deepest dip. 
Having tried each, and found no clear indication of the 
known ~4 day cycle, apart from in the case of the phase 
measure, we settled on this last method. The phase and 
error were determined as for the narrowest dips previ- 
ously. Of course in some cases the dip morphology is 
very complex, and we had to subjectively decide which 
part to fit with the Gaussian, although the generous 
FWHM error estimates will largely take this source of 
uncertainty into account. In the case of an asymmet- 
ric dip, the fit was only done on the portion which was 
roughly Gaussian in shape. The resultant dip phases 
as a function of time (the midpoint of each observing 
block) are plotted in Fig.0. 

Owing to the restricted sampling, a number of sinu- 
soids with different periods can be fitted to the data. We 
used x 2 fitting to take the dip phase uncertainties into 
account, finding local minima in the expected range at 
5.08d(x?=0.97, 10 d.o.f.), 4.74 (0.59), 4.39 (0.01), 4.10 
(0.98), and 3.84 (0.88)d. We do note that the shortest 
period fit is 3.84 ± 0.06 d which is therefore within la 
of P beat =3.90d. 

Turning to the "beat" cycle phasing of the folded 
lightcurves we still find that none of the candidate peri- 
ods yields a rcpcatable cyclic variation in morphology. 
Only for the 5.08 and 4.74 d periods, are there not in- 
consistencies within one "beat" cycle, or between the 



consecutive July cycles, as shown in the left-hand panels 
of figures H and |^. For these two, one can then at least 
achieve consistency by assuming that the "beat" phase 
can change over the course of a few cycles, as shown in 
the right hand panels. In fact the assumption of such a 
quasi-periodic "beat" period would be consistent with 
both the clock-noise known in the well-studied preces- 
sion periods of SS433 and Her X-l (Baykal, Anderson & 
Margon 1993) and the evolving supcrhump periods in- 
ferred in SU UMa systems during their super-outbursts 
(Warner 1995a). However, the shifts here are relatively 
large: for the 5.08 d period each month's data requires 
shifting relative to the others by between 0.1 and 0.5 in 
phase, whilst for the 4.74 d period the shifts are about 
0.4 in phase, but only the extreme June and September 
datasets are discrepant with July/ August. 

As a further test and to enable later comparison, 
we have also plotted the optical lightcurves (folded on 
Px) against "beat" cycle phase in figure [l0[ The phase 
has been set to agree with the (adjusted) phasing of the 
X-ray lightcurves from June. However, we see that the 
optical data does not not permit any further discrimi- 
nation between candidates. 

To summarise, we have found that Pf, ea t=3.90d is 
inconsistent with a systematic, cyclic variation in mor- 
phology. We have attempted to identify other candidate 
periods from an analysis of the varying dip phases, but 
again discrepancies are found. Only if we then relax the 
condition of strict periodicity and allow phase shifts be- 
tween well-separated cycles do we find that a period of 
4.74 ± 0.05 d can best represent the period on which the 
morphology changes repeat. 



5 LONGER TERM VARIATIONS IN FLUX 
LEVELS 

To search for any periodicities in the RXTE /ASM 
data, we calculated both Lomb-Scargle (LS) and PDM 
periodograms from the one-day average data. These 
two methods are complementary, the former modified 
Fourier transform being ideal for finding low ampli- 
tude, smooth/sinusoidal variations, and the PDM, any 
highly non-sinusoidal periodic modulations. The results 
are presented in Fig.|ll|, where the periods of interest 
are marked by dotted lines/shading, namely the previ- 
ous 199 d detection, the P^ eoi =3. 8-5.1 d range and the 
3.90 d "beat" of the X-ray and optical periods. 

However, in order to say anything meaningful re- 
garding the detection of any variability at these or any 
other frequencies we must consider the significance lev- 
els of the peaks (most reliably done for the LS peri- 
odogram). In figure ^ the log- log binned power spec- 
trum is shown, which clearly cannot be described by 
white noise, as the LS power is not approximately con- 
stant with frequency. Rather it has the appearance of a 
red noise power law, but with a break at about 0.02 d _1 
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Figure 6: Each of the X-ray datasets folded on Px and normalised by dividing by the non-dip flux level, and then 
plotted as a function of phase in the 3.90 d "beat" cycle. The folds corresponding to different cycles have been 
marked by different line styles (see key at bottom). The '?' indicate cases where the phasing is inconsistent with 
the expected systematic morphology evolution over a "beat" cycle. Right: in particular, we note the discrepancy 
between the consecutive July cycles, as replotted here for clarity. In one cycle, there are narrow dip morphologies at 
(f'beat = 0.75 and 0.98, whilst in the next cycle the same region of phase is occupied by an example of a more complex 
broad, primary + secondary dip shape (at <f>beat = 0.88). There is also a departure from systematic evolution at 
around 4>beat =0.1. 



and decreasing slope at the highest frequencies due to 
Poisson/instrumcntal white noise contributions. 

As a consequence, the method outlined in Scar- 
gle (1982) (and implemented in for example Homer, 



Charles & O'Donoghue 1998) for determining signifi- 
cance levels in the presence of white noise cannot be 
applied directly in this instance. Instead, we have de- 
veloped a modified approach, whereby an input model 
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Figure 7: Considering all the RXTE/PCA X-ray data which shows a primary dip, the phase of dip minimum is plotted 
against mid-time of observation. The periodic wandering in phase is most apparent. The candidate P{,eat = ^-7^d i s 
shown fitted to the data, along with the residuals. This fit yields an ephemeris of Tb(HJD) = 2451030.93 ± 0.27 + 
n* (4.74 ±0.05). 



noise spectrum replaces the straightforward white noise. 
An average over many resulting simulated noise power 
spectra is shown over-plotted on the data in Fig.[l2|. As 
the noise level clearly varies with frequency, we then 
split the full frequency range into 60 overlapping 0.01 
d _1 bins, for each of which we note the peak power due 
to noise in each simulated power spectrum. In this way 
we can build up the cumulative distribution function for 
the peak powers due to noise, and finally set a 95% con- 
fidence level, above which any peak in the data has less 
than a 5% chance of being due to a noise component, 
i.e the strict periodicities which we seek. 

Hence, we can see that no peaks significant above 
noise are present anywhere close to the a priori pe- 
riods of interest. Setting upper limits we find that a 
20% (semi-amplitude) at 199 d, and 22% at -4 d signal 
would be detected with less than 95% confidence. 

However, from the LS pcriodogram moderately sig- 
nificant (> 95% confidence) periodicities are easily iden- 
tifiable at about 0.012 d _1 corresponding to 83.0±0.4d 
and at — 0.043 d _1 (23.4 d ), and with smaller ampli- 
tudes at -0.083 and -O.llGd- 1 (-12, -8.6 d ), al- 



though only the first of these is evident in the PDM. In- 
deed, the 83.0 d period can be clearly seen in the binned 
lightcurve itself, as indicated in figure |[ Consideration 
of the times of maximum flux also imply that it is not 
strictly periodic. 

To shed further light on all the detected period- 
icities, we also considered the X-ray data from the 
complete band and in three energy bands folded and 
phase binned on each. The estimated peak-to-peak am- 
plitudes for all modulations are listed in table |^. For 
thc three shorter periods, the amplitudes vary in a non- 
systematic manner with energy, hence casting doubt 
upon their reality. The signal-to-noise is simply too poor 
for such low amplitude effects. Only in the case of the 
83 d signal do we see a significant modulation in each 
of the bands, although we note that its amplitude is 
much less than that of the supposed 199 d signal found 
in Vela 5B data and hence it would not then have been 
detectable. Notably, the amplitude is approximately the 
same in each band, indicating energy independence (as 
confirmed by the hardness ratios, see Fig.[l3|) . Both the 
smooth morphology and lack of significant energy de- 
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Figure 8: Left: each of the X-ray datasets folded on Px and normalised by dividing by the non-dip flux level, and 
then plotted as a function of phase of the 4.74 d "beat" cycle. The folds corresponding to different cycles have 
been marked by different line styles (see key at bottom). The '?' indicate cases where the phasing is inconsistent 
with the expected systematic morphology evolution over a "beat" cycle. Right: in order to achieve consistency, the 
datasets from June and September have been offset in phase by hand. The original ('-') and new ('+') phases for 
each discrepant dataset have been marked. 



pendence, argue against an obscuration effect such as 
that seen due to the precessing disc in Her X-l (Gerend 
& Boynton 1976; Crosa & Boynton 1980). It seems 
more likely that XB 1916-053 is similar to X1820-303, 
in which the X-ray flux shows a 176d periodicity, also 
correlated with bursting activity, and hence an inher- 



ent physical property of the system, most likely due to 
a variation in the accretion rate (see also Bloser et al. 
2000). 
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Figure 9: See as for figure g, but using the 5.08 d candidate "beat" period, and in addition the August datasets have 
also been adjusted in the right-hand panel. 



6 OPTICAL PERIODICITIES 



The optical data exhibits significant (0.5 mag) vari- 
ations from night to night over the five consecutive 
runs (see Fig. |2| lower panel) . Although larger in ampli- 
tude, they are consistent with the changes in the X-ray 
flux, as expected. In any case, for consideration of the 
~3000 s timescales, of which there are at least 2.5 cy- 
cles per night, we can safely detrend each night's data 



with a linear polynomial fit, removing any longer term 
effects. 

Once again, both an LS and PDM periodogram 
were calculated from this detrended lightcurvc (see 
Fig0). There are a number of significant peaks appar- 
ent in both periodograms. Of course the one-day alias- 
ing due to the sampling can cause confusion, but with 
a priori knowledge of which frequencies to consider, we 
may identify the principal periodicities quite easily. The 
largest peak [y = 28.570 d -1 indicated along with its 
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Figure 10: The optical data on XB 1916-053 normalized, folded on Px and binned, and then plotted at the appropriate 
"beat" cycle phasing, using the 3.90, 4.74 and 5.08 d periods respectively. The ephemcrides for each panel have been 
set so that direct comparison of X-ray and optical can be made. In the latter two cases, this phasing agrees with the 
adjusted versions. 



one day aliases) corresponds to P = 3024. 2±0. 4 s, close, 
but not consistent with the 3027.6 s optical periodicity 
of Callanan, Grindlay & Cool (1995, henceforth P opt ). 
The apparent discrepancy of 8cr is of note; it may sim- 
ply indicate the underestimation of the true uncertainty 
in the peak frequency, due perhaps to the effects of the 
larger amplitude (and varying) modulation present on 
the nearby X-ray period, or indeed may be evidence 
for a period change. We will consider the stability of 
the optical period further in the next section. Exclud- 
ing these, the next most significant peak is indeed at 
2995.8 ± 0.6 s (28.84 d" 1 ), which is close to the X-ray 
period. Subtraction of the longer period modulation re- 
veals the signal on Px even more clearly. Notably, this 
has a revised period of 3000.2 ± 0.6 s consistent with 
our X-ray determination. Moreover, this signal appears 
relatively more distinct in the PDMs, implying that it 
is non-sinusoidal (as one would expect). Such an occur- 
rence of both periodicities in the optical was also seen 
in data from 1990 September by Callanan, Grindlay & 
Cool (1995), when they comment that XB 1916-053 was 



in an anomalously low optical state, and by Grindlay 
(1989). 

As we did for the X-ray data, we then proceeded 
to examine the folded and binned lightcurves for each 
day and the complete dataset, using both the period 
and ephemeris from Callanan, Grindlay & Cool (1995) 
and our own X-ray period. The overall lightcurve folded 
on P op t shows an almost sinusoidal profile, but with a 
much smaller amplitude (~10%), than the nightly folds 
(up to ~25%). In common with the X-ray folds, there is 
clear evolution of the morphology from day-to-day. In 
particular the phase of the principal minimum appears 
to shift systematically relative to the plotted sinusoid. 
If, indeed, the 3027s period represents that of the binary 
orbit, so long as the reprocessing geometry is in some 
way asymmetric, the change in its aspect will lead to 
the smooth modulation we see. Hence, we fitted and 
subtracted a sinusoid fit to the overall fold (as a good 
approximation to its shape), to leave the effects presum- 
ably due to obscuring disc structure alone. The folds in 
Fig. [l^ shed further light. It is now clear that the re- 
maining structure in the optical lightcurve varies on Px 
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Figure 11: Phase Dispersion (upper) and Lomb-Scarglc (lower) periodograms of the complete ~4 yr RXTE/ ASM 
lightcurve of XB 1916-053. A number of peaks are present, but there are none corresponding to either the previously 
detected 199 d modulation or the expected 3.90 "beat" period (marked by dashed vertical lines) nor any even in the 
range of possible "beat" periods (shaded region). However, those at ~0.012, ~0.043, ^0.083 and ^0.116d _1 (~83, 
~23, ~12, ^8.6 d ) are statistically significant (>95% confidence) in excess of the red-noise level in their region of 
the power spectrum. 



not P op t This strongly infers that the same effect caus- 
ing the X-ray dips, creates corresponding optical 'dips'. 
Comparing X-ray and optical folds (on Px) from the 
same phase of the "beat" cycle (Figs. p|-|lO|), we see that 
there is surprisingly little evidence of a clear relation- 
ship between the two. For instance, in the simultaneous 
data of June 25 the X-ray exhibits a distinct narrow 
dip, whilst there is very little optical modulation at all. 
Nevertheless, we may determine the phase relationship 
of the dips in each band in the other cases, finding that 
the optical dip lags the X-ray by ~0.2-0.3 in phase. 



6.1 Revised optical ephemeris 

We may also reconsider the optical ephemeris derived 
by Callanan, Grindlay & Cool (1995), by using the 
complete optical lightcurve folded on their ephemeris 
(see Fig. [l5|a lowest folded lightcurve) to extend the 



baseline by a further 4 years. Clearly, the minimum is 
shifted relative to phase zero. Applying their method 
we fitted a Gaussian to the principal minimum giving 
4>min — — 0.20±0.09. The a priori probability of obtain- 
ing a minimum this close to the published ephemeris 
(equivalent to ~ler) is still ^40%. Hence, the period 
could have changed significantly since the earlier obser- 
vations, as suggested by the periodogram result, and 
we simply have a chance approximate alignment. On 
the other hand, if we assume that the period is stable 
then the shift implies that the 1995 ephemeris is slightly 
incorrect. So to update it, we also digitised the residual 
values as presented in Callanan, Grindlay & Cool (1995) 
Fig. 6 (but with the corresponding HJD values calcu- 
lated from the quoted observation time/dates). In fig- 
ure |l^ we present all the residuals to P op t = 3027.551 s 
from 1987 through 1998. If this period is correct then 
a constant fit will be acceptable, whereas a slightly 
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Figure 12: Log-Log binned LS periodogram (fine histogram, with error-bars) of complete RXTE/ ASM lightcurve 
of XB 1916-053. Note the general form, with approximately a power law between the break at ~0.02d _1 and 
0.1 d _1 (red noise part), which tails off to higher frequencies as the white noise contribution (flat spectrum) becomes 
dominant. The components of the input model noise spectrum (white noise, plus red noise, with a break at 0.02 d _1 
) are shown offset vertically and the resulting average (over 450 simulations) is over-plotted (in a thick line). 



Table 2: Semi-amplitude of flux modulations in the RXTE / ASM lightcurves 



Energy range 


Mean ASM 


Percentage peak-peak flux 


variation 


(keV) 


Count rate (s _1 ) 


83.0d 23.4d 12. Id 


8.6d 


1.3-3.0 


0.44 


60 ±20 10 ± 10 30 ±20 


50 ±20 


3.0-5.0 


0.38 


50 ± 20 20 ±20 10 ±10 


10 ±10 


5.0-12.1 


0.54 


50 ± 20 50 ± 20 20 ± 20 


20 ± 10 


1.3-12.1 


1.16 


60 ± 10 40 ± 10 20 ±10 


20 ± 10 



different period will be represented by a linear term. 
For the constant fit we find an acceptable xt = 1-06, 
but with a revised T (HJD) = 244, 6900.0019 ± 0.001, 
whilst for a new period and ephemeris of Tb(HJD) = 
2446900.0128 ± 0.0055 + n * (0.03504115 ± 0.00000002) 
{P' opt = 3027.555 s) we have a smaller X l = 0.77. 
Applying a one-sided F-test (Bevington & Robinson 
1992) confirms that the extra linear term does provide 
a marginally better fit than the constant with 91% con- 
fidence. However, the revised period is only la larger 
than the Callanan, Grindlay & Cool (1995) result, hence 
not too much weight should be attached to this mi- 
nor revision. The most important aspect of this new 
ephemeris is that our data, on face value, is consistent 
with all previous observations, that is with the stabil- 
ity of the optical modulation over more than a decade. 
If confirmed this would provide yet stronger support 
for its orbital origin, but still further observations are 



needed to do so, preferably spaced by a few months to 
maximise the accuracy of the cycle count. 

6.2 Further comparison of optical and X-ray 
modulations 

Having considered the optical lightcurve folded on the 
X-ray period, we will now look at the converse. In fig- 
ure [I?] we present the complete X-ray lightcurves folded 
(and phase binned) on both Px and P op t- The presence 
of both primary and secondary dips becomes clear in 
the Px fold, showing how the latter can be as deep, 
but only occur a fraction of the time. Before folding on 
Popt we subtracted the mean modulation on Px (i.e. 
the binned fold) in the same way we treated the opti- 
cal data. When folding on this slightly longer period, 
the primary dips are spread over a much larger phase 
region, and are only limited due to the effects of sam- 
pling; leading to the apparent modulation in the phase 
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Figure 13: Left: One-day averaged RXTE/ ASM data on XB 1916-053 in three energy bands and complete folded and 
phase binned on an 83 d period. The lightcurves have been normalised and then offset vertically (by 1 cts s -1 ) for 
clarity. Right: hardness ratios between the different bands clearly demonstrate the approximate energy independence 
of the modulation amplitude at this period (right). 



binned plot. Lastly, considering only the non-dip data 
we see that there is in fact no evidence of any overall 
modulation unrelated to the dip behaviour. 



7 DISCUSSION 

Taking a qualitative approach we will consider the im- 
plications of our results. A range of models have pre- 
viously been put forward to explain the observational 
properties of XB 1916-053 ranging from a hierarchical 
triple system, to something akin to the Cataclysmic 
Variable (CV) superhumping systems with their pro- 
cessing discs. We also outline a modification to the lat- 
ter, related to the more recent discovery of negative 
supcrhumps. Dealing with each model in turn, we will 
confront its predictions with our detailed observational 
results. 



7.1 Hierarchical triple 

Grindlay et al. (1988) and Grindlay (1989) proposed 
that XB 1916-053 is part of a hierarchical triple system 
in order to provide an explanation for the two periodic- 
ities close to 3000 s and the then reported longer-term 
199 d flux variation. In their model, the ultra-compact 
3000 s binary, whose orbital period is that measured 
in the optical, is orbited by a third star on a some- 
what longer, retrograde orbit (P or b — 4d). The gravita- 
tional influence of this third star leads to perturbations 
of the mass transfer rate on the X-ray period. Unlike 
other dippers, they require that the dips arise when 
the mass transfer rate is periodically enhanced, thereby 
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Figure 14: Phase Dispersion (upper) and Lomb-Scarglc (middle) pcriodograms of the complete 5 night optical 
lightcurve of XB 1916-053 (left) and the same data after subtraction of the measured sinusoidal "orbital" modulation 
at 3024.2 s (right). The 1 day aliasing due to the sampling complicates the periodograms, but as shown with the 
simulated data (bottom) the structure is reproducible simply by supposing there are two signals at close to Px 
and P op t- Note: (left) the measured peaks are shifted relative to the known Px = 3000.6s and P op t = 3027.6s as 
indicated by the vertical lines (dashed and dotted), whilst after the removal of the latter (right) the measured Px 
agrees with our X-ray determination. 



forming extended disc-edge structure which then cir- 
cularises (around the disc) within one binary orbit. In 
this way, the problem of the dip-causing structure be- 
ing fixed in binary phase is solved. Moreover, the long- 
term flux modulation is also explained, as arising from 
the slow changes in binary eccentricity at the period 
Piong ~ Pouter/ Pinner (Mazeh & Shaham 1979), where 
Pinner refers to the inner binary orbit and P ou ter to the 
outer retrograde orbit of the third star. 

However, one must bear in mind, before even con- 
sidering the observational constraints, that the forma- 
tion of such a hierarchical triple (with a retrograde third 
star) is only probable within the crowded environment 
of a globular cluster core. Indeed Grindlay (1989) sug- 
gests that this may have been the case for XB 1916-053. 

7.1.1 Comparison with observations 
In favour: 



(i) Popt — Porb hi agreement with the apparent high 
stability of the overall optical modulation. 

Against: 

(i) There is no reason to suppose XB 1916-053 has a 
different dip formation mechanism to the other dippers 
and hence we might expect some evidence for X-ray dips 
at P or b, in addition to the dips occurring at times when 
the mass transfer is enhanced (unless an extremely spe- 
cific orbital inclination is chosen so that orbital effects 
are never visible). 

(ii) There is no confirmation of the old 199 d period, 
and the new 83 d value is possibly too short to fit in with 
changes in eccentricity driven by the third star. Follow- 
ing Grindlay (1988) we have P [ong = K x P^ uter / P inner , 
where now Pi ong = 83 d, P !lmer = 0.035 d and P ou ter ~ 
| x P beat - | x 3.90 - 2.6 d, yielding K = 0.4 com- 
pared to their 0.7 for a constant which is expected to 
be of order unity. 
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Figure 15: a) The optical lightcurves from each night and overall folded on P op t (grey points) and phase binned 
(over-plotted in black). The data was firstly normalised by dividing by the mean instrumental flux value, and they 
have been plotted vertically offset for clarity. Note that the nightly modulations are varied in morphology and share 
no phase relation to each other or to the overall fold, which is also much more smoothly sinusoidal in appearance, 
b) Same as (a) except that the 3000.6s X-ray period has been used for folding. Now the clearest features of the 
modulation appear at the same phase, indicating that this shorter period is the dominant, over the course of a few 
cycles at least (c) Before folding on Px, the sinusoidal fit at P opt (good approximation) to the complete optical 
dataset has been subtracted. The features which must be related to the X-ray dips are still clearly present. 



7.2 Permanent superhumper 

The presence of doubly-periodic optical modulations is 
a well-known phenomenon within the field of CVs. The 
effect was first identified during the super-outburst 
events of the SU UMa sub-class of CVs, and hence 
the non-orbital modulations were referred to as su- 
perhumps (Vogt f984; Warner 1975). Further study 
of CV lightcurves led to the discovery of permanent 
superhumping systems, amongst the nova-like and 
nova remnant classes (Warner 1995b). In all cases, 
the systems have extreme mass ratios, which in turn 
leads to the tidal distortion of the disc by the donor 
(Whitehurst 1988; Whitehurst & King 1991). The 
elliptical disc then precesses in a progradc manner 
(in the inertial frame) on a timescale of days. Both 
White (1989) and Smalc et al. (1992) have proposed 
that XB 1916-053 is an example of an analogue LMXB 
system, although the authors present differing ideas as 
to the origin of the optical modulation. White (1989) 
invoked a purely geometrical effect, whereby the pro- 
jected area of the irradiated disc (the principal optical 
emitter in the system) varies as the disc precesses. 
However, the successful Whitehurst (1988) model for 



supcrhumps shows that such a variation would be on 
the precession period, not on its beat with the orbital. 
Smale et al. (1992) presented a modified approach, in 
which the vertical extent of the pressure-supported disc 
would vary on P SU perhump, in response to the varying 
vertical component of gravity due to the secondary. 
Hence, in turn the area available to reprocess X-ray 
emission into optical would also vary, leading to the 
optical modulation we observe (also see O'Donoghue & 
Charles 1996). Indeed, it has been proposed that this is 
the analogous mechanism that works for CVs, except 
that it is the increased reprocessing of ultraviolet 
emission from the innermost disc which provides the 
superhump light source (Harlaftis et al. 1992, see also 
Billington et al. 1996). Previously, the direct enhance- 
ment of viscous energy release in the outer disc had 
been assumed to provide the light source (Whitehurst 
1988). In any case, such superhump periods will always 
be longer than the binary period (by a few percent). 
Hence, in contrast to the triple model it requires the 
X-ray period to be that of the binary and the optical 
due to the superhump effect. 
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Figure 16: Using the Callanan, Grindlay & Cool (1995) optical period the residuals between the predicted time of 
optical minimum and that measured is plotted for all the observations 1987 through 1998, together with a constant 
fit, as appropriate for their ephemeris. Clearly, our 1998 point has quite a large residual implying the possible need 
for a longer period (as indicated by the linear fit). Nonetheless, note how small the scatter appears over the entire 
11 year baseline, inferring its high stability. 



7.2.1 Comparison with observations 
In favour: 

(i) With Px = Porb the X-ray obscuring material 
is straightforwardly fixed within the binary frame, and 
the same mechanism can work as for other X-ray dip- 
pers. The dominant optical modulation on this period 
is also readily explained. Both the vertical structure of 
the impact region and the remainder of the disc will 
be emitting optically due to reprocessing of intercepted 
X-rays, and hence obscuration effects and changes in 
the aspect of the disc structure can account for that 
modulation. 

(ii) Presuming this material to be related to the 
stream impact region at the disc edge, the presence of 
an elliptical disc provides a ready explanation for the 
smooth variation in primary dip phase over the preces- 
sion cycle. 

(iii) From previous observations, the longer period 
optical modulation dominates over the variations on Px 



when the system is brighter, and indeed this is exactly 
what we would expect if the superhump light source 
is due to irradiation (in contrast to the comment of 
Callanan, Grindlay & Cool 1995) 

Against: 

(i) The very distinct dip morphology changes, and 
especially the disappearance of primary dips altogether 
is difficult to explain simply with this model, but this is 
probably not a major drawback as more detailed mod- 
elling is required. 

(ii) Ongoing work on the permanent CV supcr- 
humpers does imply that even in these system with ap- 
parently stable mass transfer and hence disc size (c.f. 
the outbursting SU UMas) the superhump periods are 
almost stable over the course of weeks, but not years. 
This poses a problem with interpreting the very sta- 
ble optical modulation as some sort of superhump, i.e. 
with putting P opt = P SU perhump, given the very long 
timescale of apparent stability of XB 1916-053. 
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Figure 17: Upper panel: Complete June-September RXTE/PCA lightcurvc on XB 1916-053. Prior to folding, each 
of the daily datasets was normalised by dividing by the non-dip flux level, a) Lightcurve folded and phase binned 
on the X-ray period of 3000 s. b) After subtraction of this average (binned) modulation from the lightcurve, folded 
once more on Px (offset by 1 for clarity), c) Same subtracted data as in (b), but now folded and phase binned on 
the optical period of 3027s (offset by 1 for clarity). An apparent modulation on P op t remains even after subtraction 
of the mean X-ray modulation shown in (a). Lower panel: Only the non-dip data (i.e. applying a 75% flux cut to 
remove the dip data points). Without the dip data there is no longer any modulation on P op t, showing that this only 
arises from the uneven sampling of the dips. 



7.3 New model: permanent negative 
superhumper 

In addition to the CV superhumps with periods longer 
than orbital, (positive), there is now evidence for such 
modulations with periods a few percent shorter, known 
as negative superhumps (Patterson et al. 1993). These 
often appear together with the positive kind, and in 
fact the proposed explanation is essentially the same, 
except that one now requires retrograde precession of 
some feature associated with the disc (e.g. the line of 
nodes). However, this will only in general arise if the 
disc can be made to lie outside of the orbital plane, i.e. 
it is tilted and/or warped, which is the main difficulty. 
In the case of LMXBs, we do at least have some idea 
of how to form such stable warped discs. Irradiation- 
driven coronal winds have been put forward by Schandl 
k, Meyer (1994) in the case of Her X-l, where as soon 
as the disc is perturbed out of the plane, the torque on 



material leaving the disc can act to enhance and then 
preserve the warp. We note that XB 1916-053 does in- 
deed possess an accretion disc corona, supporting the 
presence of such a mechanism. Another possibility could 
be magnetic effects, as posed by Lai (1999), although 
this is only likely to be of importance in more highly 
magnetic systems and in any case close to the magnc- 
tosphcrc, i.e. the innermost disc regions. Unfortunately, 
further work by Ogilvie & Dubus (2000) on the appli- 
cation of irradiation-driven warping to X-ray binaries 
has convincingly shown that this will not work for most 
LMXBs and certainly not for such an ultra-compact sys- 
tem as XB 1916-053. Nevertheless, as in the case of CVs 
let us postulate that XB 1916-053 possesses a warped, 
retrograde, precessing disc, and see what this model will 
do for us. 

Considering the optical emission, the underlying 
modulation (on P op t) will arise from some asymmetry 
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in the reprocessing geometry, which then leads to a con- 
tribution which varies with binary phase as the aspect 
changes, for instance the heated face of the secondary. 
However, XB 1916-053 has both a very small compan- 
ion and significant vertical disc structure, hence, as was 
found in our modelling of the multi-colour lightcurves of 
4U1957+115 (Hakala, Muhli & Dubus 1999), shadowing 
could easily render its modulated contribution negligi- 
ble. This same shadowing could be partly responsible 
for the dominant shorter (X-ray) period optical modu- 
lation, together with both enhanced emission from and 
self-obscuration of the disc, all related to the negative 
superhump effect itself. 

Moving on to X-ray dips: as in the standard for- 
mation scheme, they arise from the obscuration of the 
central X-ray source by material thrown up at the im- 
pact point of the accretion stream and the disc. But 
what is novel is that the radial and azimuthal position 
of this region varies systematically with the precession 
of the disc. This was first put forward by Schandl (1996) 
to explain the pre-eclipse dips seen in Her X-l, which 
arc definitively known to recur on a period less than the 
orbital (Scott & Leahy 1999). The warped disc, with a 
retrograde, preccssing line of nodes, allows the accre- 
tion stream to overshoot the disc edge itself for a large 
fraction (but not all) of the precession cycle, leading to 
dips which occur progressively earlier in binary phase. 



7.3.1 Comparison with observations 
In favour: 

(i) With P opt = P or b as was the case in the triple 
model, the stability of the optical modulation is consis- 
tent with its orbital origin. 

(ii) There is the possibility to account for the sys- 
tematic dip morphology changes. Due to the warping 
the accretion stream will impact the disc at varying an- 
gles of incidence and depths within the neutron star's 
potential well (affecting the energetics). 

Against: 

(i) As a corollary to the warp leading to dip shape 
changes, we encounter the problem that as soon as the 
disc leaves the mid plane, dips can only occur for half 
the precession cycle when the stream intersects the face 
of the disc we see. So one can explain the absence of 
primary dips easily enough, but for far too much of the 
cycle (as in Her X-l). 



7.4 Summary 

Following the comparison of our observational data with 
a number of possible models for XB 1916-053, we find 
that both models invoking the optical period as that 
of the binary have serious problems. In the cases of 
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the hierarchical triple model and the permanent neg- 
ative superhumper inconsistencies are found with their 
predictions for X-ray dip formation. The most convinc- 
ing model is that of a permanent positive superhumper, 
where a tidally distorted, elliptical disc precesses in a 
prograde manner, leading to the two close periodicities, 
and requiring the X-ray dips to recur on the binary 
orbital period. However, this promising model can not 
explain all the observational features. Specifically, more 
work is needed to confirm whether indeed the X-ray dip 
period is more stable than the optical (consistent with 
its proposed orbital origin) . This will require reanalysis 
of archival X-ray data, and further observational effort 
optically. Lastly, this model provides no ready explana- 
tion for the dip morphology cycle and indeed none of the 
models can explain the discrepancy between its 4.74 d 
period and the 3.90 d "beat" period of ~3000s signals. 
Clearly, detailed modelling of the combined correspond- 
ing X-ray and optical orbital modulations is needed to 
tackle these issues, but such work is beyond the scope 
of the current paper, and will be presented in a forth- 
coming one. 
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